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We offer predictions of symmetron modified gravity in tlie neighborhood of realistic dark matter 
halos. The predictions for the fifth force are obtained by solving the nonlinear symmetron equation 
of motion in the spherical NFW approximation. In addition, we compare the three major known 
screening mechanisms: Vainshtein, Chameleon, and Symmetron around such dark matter sources, 
emphasizing the significant differences between them and highfighting observational tests which 
exploit these differences. Finally, we demonstrate the host halo environmental screening effect 
("blanket screening") on smaller satellite halos by solving for the modified forces around a density 
profile which is the sum of satellite and approximate host components. 
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I. INTRODUCTION 

The observed acceleration in the expansion of the uni- 
verse can arise from a dark energy component or from 
a departure of gravity from general relativity (GR) on 
cosmological scales. One way to distinguish between the 
two possibilities is to consider the growth of perturba- 
tions. For modified gravity (MG) theories, the relation 
of the expansion history to the growth of perturbations 
is specific to every model. In the quasi-static, Newto- 
nian linear regime, several authors have parameterized 
the growth of perturbations with g{k,z) = G/Gn and 
il{k,z) = i;/<P (e.g. m\)- 

Laboratory and solar system constraints (see [4j and 
referrences therein) require any viable MG theory to 
have some mechanism by which it mimics the predic- 
tions of GR within Milky Way-size halos. Such "screen- 
ing" mechanisms (see [5] for a review) generally deter- 
mine the deviations from GR based on the local den- 
sity: in high density environments the scalar force is sup- 
pressed, while in low density environments it can be of 
approximately gravitational strength. Two such mecha- 
nisms have been extensively studied in the literature: the 
chameleon screening [B^ of f{R) theories [7] and Vain- 
shtein screening [8^ of higher dimensional (e.g., Dvali- 
Gabadadze-Porrati (DGP) gravity [S]) and Galileon mod- 
els 10]. Recent work has included detailed simulations 
which arc necessary because of the nonlinearity inherent 
in how GR is recovered inside the Milky Way [TTMlSj . 

In this work we explore the symmetron model of |14j 
(see also [T5] and [16]), which exploits a novel screening 
mechanism similar in part to chameleon screening but 
with key differences. The symmetron has a vacuum ex- 
pectation value (VEV) that is large in low density envi- 
ronments and small in high density environments. Sym- 
metron screening then relies on a coupling to matter that 
is proportional to the VEV, thus suppressing the scalar 
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force in high density environments. Recent work has 
focused on symmetron cosmology, including the evolu- 
tion of the symmetron field through various cosmological 
epochs [T7] , as well as its effect on linear (TB] and nonlin- 
ear |19j structure formation. 

Tests of gravity on linear scales have some limitations. 
The g, rj parameterization is only valid on scales smaller 
than the superhorizon regime and larger than the nonlin- 
ear regime. At high redshift there is a different problem: 
since MG models recover GR at high redshift for consis- 
tency with CMB and Nucleosynthesis observations, ef- 
fects of enhanced forces are manifested only at late times. 
Thus, even if observations are made late enough to be 
within the MG era, the signal has had limited time to 
accumulate. 

In this study we consider tests on scales within and out- 
side virial radii of dark matter halos modeled with the 
Navarro-Frenk- White (NFW) [2^ profile. In this regime 
the predicted deviations due modified gravity can be sig- 
nificantly larger than in the linear regime (measurement 
errors and systematic uncertainties need to be taken into 
account but will not be considered here). As highlighted 
by [3Tj and [55^, astrophysical tests in this regime can 
provide effective tests of chameleon theories. We will 
calculate the predicted deviations for symmetron theo- 
ries. 

In § II we describe the symmetron theory and our 
method for calculating modified forces around NFW ha- 
los. § III contains our isolated halo results as well as a 
comparison of the various screening mechanisms. In § IV 
we model and show results for forces on test particles in 
two-body host-satellite systems. We conclude in §V. 



II. FORCE PROFILES OF NFW HALOS 
A. Symmetron Theory 

In the Einstein frame we can describe the gravitational 
forces as GR with an additional "fifth force," mediated 
by the symmetron field, 4>. For GR we have the usual 
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Poisson equation: V^^'n = ^T^Gp, leading to 
d^-N GM{< r) 



dr 
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The symmetron equation of motion in the presence of 
non-relativistic matter |14| is 
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Note that the relevant field range is <C Afg, such that 
any 0{(j)'^/M^) terms in A{(j)) can be consistently ne- 
glected. The potential V{(j)) comprises the most general 
renormalizable form invariant under the Z2 symmetry 
(j) — > —(p. The coupling to matter ^ (fp' /M"^ is the leading 
such coupling compatible with the symmetry. The model 
involves two mass scales, ^ and Mg, and one positive di- 
mensionless coupling A. The mass term is tachyonic, so 
that the Z2 symmetry (/)—>■—(/) is spontaneously broken. 
The effective potential of Eq. ^ is 
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Whether the quadratic term is negative or not, and hence 
whether the Z2 symmetry is spontaneously broken or not, 
depends on the local matter density. 

The screening mechanism works roughly as follows: in 
vacuum or in large voids, where p ~ 0, the potential 
breaks reflection symmetry spontaneously, and the scalar 
acquires a VEV \4>\ — 4>o = m/V^; in regions of high 
density, such that p > M^fi^, the effective potential no 
longer breaks the symmetry, and the VEV goes to zero. 
Meanwhile, to lowest order the symmetron-matter cou- 
pling is ^ pcfy^/M^. Fluctuations Scj) around the local 
background value (j>vEV, which would be detected by lo- 
cal experiments, couple to density as 
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In particular, the coupling is proportional to the local 
VEV. In high-density environments where the symmetry 
is restored, the VEV should be near zero and fluctuations 
of (j) do not couple to matter. In less dense environments, 
where p < M^p'^ and the symmetry is broken, the cou- 
pling turns on. 

For a static-spherically symmetric source, Eq. ^ be- 
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We set the parameters as in [T3]: Mg = 10~^Mpi sat- 
isfies solar system constraints while still allowing for or- 
der unity deviations elsewhere. Also, p = y^/Mg and 
A = (/i/0o)^j where Pc is the average cosmological density 
today and (po is the background value of the field. Note 
that these parameter choices correspond to p ^ Mpc~^, 
constraining symmetron effects to ~ Mpc distances. 

As in [17] . the symmetron- mediated force relative 
to the Newtonian force between two test masses in 
vacuum is set by the symmetry-breaking value 00 : 
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If the scalar-mediated force is to be comparable to gravity 
in vacuum, then we must impose (po/AI^ ^ 1/Mp\, that 
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where 5 ^ 0(1). To be precise, it follows from Eq. ([8]) 
that g measures the strength of the scalar force in vac- 
uum relative to gravity: — 2g^Fj^. For comparison 
to f{R) and DGP theories, for which the fifth force is at 
most 1/3 i^Ni we will set g = l/Vd; otherwise we choose 
5 = 1, which is still consistent with solar system tests, 
as shown in [T3j. Note that Eq. ([s]) has no dependence 
on the test bodies involved. Extended mass distributions 
affect the scalar and Newtonian forces differently so that 
solving for the Newtonian potential and scalar field pro- 
file is required in order to evaluate the ratio. The ratio 
of forces on a test mass in the neighborhood of such an 
extended distribution is 



F^ ^ {4>/AQ{V^/M,) 
Fk V*n 



B. NFW Halos 
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We consider gravitational forces in the neighborhood 
of NFW halos ,201 , whose density profiles are a good fit 
to those of stacked simulated halos. The (untruncated) 
density is 
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where Ps and are parameters that depend on halo mass 
(see below). We define the mass M300 of each halo as 
that enclosed within the virial radius, i?3oo- This is the 
radius at which the average density enclosed is 300 times 
the critical density p^. The concentration connects the 
virial radius to the scale radius, c = R^Q^/rs- We take it 
to be a function of halo mass. 
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in the case of a top-hat density profile 
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FIG. 1: Force deviation F^/F-f^ for a top-hat density pro- 
file of total mass lO^^M© calculated numerically (dashed 
line) and analytically (dotted line). We find similar agree- 
ment between the two methods throughout the mass range 



as found in |23| . The mass and concentration definitions 
are chosen to allow comparison of our results with the 
chameleon model of [23]. We note that the profile itself 
will change somewhat with the modified forces. 
The mass enclosed within radius r is 



M(< r) = Ma, 
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where F{x) = In (1 -|- x) — x/ (1 + x). We use Eq. (13) in 
^ to solve for the Newtonian force. In order to obtain 
the symmetron profile, it is useful to define a dimen- 
sionless scalar field, ip = <j>/4'o whose equation of motion 
follows from Eq. ([t]): 
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We obtain the radial profile by substituting Eq. (Ill in 
( 14 1 with boundary conditions 
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dr 



= 0, ijj(r — )• oo) = 1 . 



(15) 
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We use a shooting algorithm to solve this nonlinear equa- 
tion, tuning the boundary condition ip(r = 0) such that 
the field stays within 1% of ip = 1 at large r for at least 
25 virial radii. Note that Eqs. (14), (15) depend only on 
the theory parameters /i and Mg] they are independent 
of 00 and A. The only effect of changing 0o is to set the 
overall amplitude of the symmetron profile and therefore 
the amplitude of the ratio of forces. 

For comparison to the work of [S] and as a check of our 
numerical solutions, we also solve for the force deviation 
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where po — 300pc- For the top- hat profile, specifying the 
mass and density fixes i?3oo- To obtain an analytic solu- 
tion, we approximate the symmetron equation of motion 
as quadratic around the appropriate minimum inside and 
outside the object (see [H] for details) resulting in 



0i„(r) = A^sinh ^r^^p^ 
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and solve for A and B by matching at the boundary. The 
exact solution is obtained using our shooting algorithm 
as in the NFW case. Figure [l] shows the equivalence 
of these methods for a top-hat mass 10^'^ M©. We find 
similar agreement between the two methods throughout 
the mass range 10^° - 2 x IQ^^Mq. 

We can define a . as an average of the force devia- 
tion over the virial radius of a halo (see |24] for details). 
This quantity can be determined from both theory and 
observations. From the theory we have calculated 
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where the integral is over the virial radius i?3oo- This is 
straightforward to compare to observations, which yield 
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where Afaoo.dyn is a dynamical mass, and A/300 is the 
"true" or lensing mass. 

As will be seen later, the deviation from GR is signif- 
icant out to ~ 10 times the virial radius. Therefore we 



define a second average £d exactly as in Eq. (18), except 
the integrals are taken over distances d = 0.5, 1, 4, and 
10 Raoo- 



III. RESULTS FOR ISOLATED HALOS 

Figure [2] shows the force deviation /Fn on a test 
particle given by Eq. ( 10 ) in the neighborhood of isolated 



NFW halos of various masses. Also pictured are the devi- 
ations for constant density spheres of the same mass. We 
plot using the symmetron vacuum value 0o — /M-p\, 
but as argued in Sec. |llj letting (/)o g(t>vi simply shifts 
these curves down by a factor . For all but the smallest 
halos (lO^^M©), the two profiles give significantly differ- 
ent results all the way out to the virial radius i?3oo- 

In the left panel of Fig. |3j we plot g^.^^ given by Eq. ( 18 ) 
for three values of the symmetron vacuum value: 0o = 
gM^/Mpi with g = 1,1/73,1/^6. Setting g = 
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FIG. 2: Ratio F^/F^ of the symmetron mediated force to the GR force for halos of different masses with field symmetry 
breaking value 0o = M^/Mpi (i.e. g = 1). The halos are modeled with NFW (solid line) and spherical tophat (dashed 
line) profiles. Note that the ratio is plotted vs. radius in units of each halo's virial radius, -R300. The background compton 
wavelength, ~ 1 Mpc (vertical dotted line) gives a sense of the physical distance. 




FIG. 3: {left panel): Averaged force deviation, ^^j^., as a function of "true" or lensing mass for different values of the symmetry 
breaking field value, (t>o = <?Mf/Mpi. {right panel): A modified averaged force deviation for the smallest field value g — l/\/6. 
The average is taken over a distance d from the center of the halo, with d — 0.5, 1, 4, and 10 -R300 from bottom. On both 
panels, the horizontal dotted line at 4/3 shows the maximum average for f{R) modified gravity. 
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FIG. 4: {left panel): Here we plot £^ = 1 + F^/Fs^ for comparison with gf(R) in iMj- Solid curves are from our symmetron 
model with 4>o ~ /VGMpi. Dashed curves are from normal branch DGP with crossover radius = 500 Mpc, = 3100 
Mpc, and self-accelerating DGP, from top to bottom, (right panel): Reproduced from [23]. Force deviation for Hu-Sawicki 
f{R) around NFW halos. Arrows denote the radius where the chameleon thin-shell condition is first met. See [24] for details 
of the DGP and f{R) models. We see that the three types of screening predict distinct transitions in -F^/Fn with respect to 
both radius and mass (see text for details). 



fixes the max deviation at 4/3 as in f{R) and DGP theo- 

ries. Equation ( 18 ) implies Mdyn,3oo — 3^-„ -^300 so that 
the symmetron theory predicts, e.^., the dynamical mass 
of a 10"Mq halo to be 50%, 25%, or 10% greater than 
the lensing mass for the three pictured values of g. 

The right panel of Fig. [3]shows the same average taken 
out to larger radii d x i?3oo, corresponding to 0.5, 1, 4, 
and 10 times the halo's virial radius. Although the peak 
deviations in Fig. [2] are at approximately 7 i?300i the 
weighting by density in the integral of Eq. ( 18 1 results in 
a relatively small increase in ^^.^ with d for a given halo. 

In Figures |4] we collect results for screening in the sym- 
metron, DGP, and chameleon theories. The DGP and 
chameleon results are from [25 ■ The left panel overlays 
our symmetron results (setting g = l/-\/6) with three 
DGP models exhibiting Vainshtein screening. These are 
normal branch DGP with crossover radius Tc = 500 Mpc, 
Tc = 3000 Mpc, and self-accelerating DGP. The right 
panel of Fig. |4] shows chameleon screening in the f{R) 
model of [7| with field cosmological value |/i?o| = 10~^. 
We refer the reader to [24 ioi further details of these 
theories. 

The "thin-shell" effect of the chameleon, in which 
only a thin shell at the edge of a screened object con- 
tributes to the fifth force, is evident in the rapid rise of 
the force deviations from zero to their peak at r « 2i?3oo. 
In contrast, partially screened halos in the symmetron 
model show nonzero deviations at smaller radii that in- 
crease all the way to ~ 7i?3oo before declining back to 



zero (see Fig. [2] for the large r behavior). DGP models 
also show nonzero deviations that increase beyond 2i?3oo, 
eventually approaching a constant value. 

The transition of the force deviations with halo 
mass also shows promise for distinguishing between the 
three screening mechanisms. For concreteness, con- 
sider the deviations at the virial radius i?3oo- For the 
chameleon background field value l/^jol — 10"^ (10~^), 
a 1O^^M0 (IO^^Mq) halo shows no deviation while a 
10^'^ M0 (IO^'^Mq) halo exhibits the maximum allowable 
deviation of 4/3. All smaller halos will likewise show a 
4/3 deviation. In contrast, the symmetron 2 x 1O^^M0 
halo is completely screened at -R300 while the deviation 
is only maximized at 4/3 for dwarfs of mass < IO^Mq. 
Thus, the degeneracy in the deviation exhibited between 
smaller halos in the chameleon model is not present in the 
symmetron. Furthermore, DGP models show no depen- 
dence on halo mass. These arguments imply that probes 
of modified forces spanning the mass range 10^ — IO^^Mq 
would be effective at distinguishing between all three 
types of screening. 



IV. HOST-SATELLITE EFFECTS 
A. Model 

We now consider the force deviation that would be ex- 
perienced by a test mass in the neighborhood of a satel- 
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FIG. 5: We average the host density in spherical shells about 
the satellite's center. The parameter d is the center-to-center 
distance between the halos, r is the standard radial coordi- 
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FIG. 6: Host-satellite density profile as a function of radius 
in units of the satellite's virial radius, centered at the satel- 
lite's location. Pictured are the NFW satellite profile (dot- 
ted), averaged host profile given by Eq. ( |20[ ) (dot-dashed), and 
the sum of host and satellite components (solid). The latter 
two curves are shown for host-satellite separation d — 1-Rhost 
(thick lines) and d = 2iihost (thin lines). In each case, the 
cusp in the total profile is located at the separation distance 
d. This plot assumes Af^at = IO^^Mq and Mhost = 2 x 10^", 
although the profile is qualitatively similar for all relevant 
satellite and host masses. The host density becomes the dom- 
inant component at r ~ Raa.t- 



lite halo which is itself blanket screened by a nearby host. 
We model host and satellite with NFW profiles as before, 
determined by Eqs. ( 11 ), ( 12 ). To preserve spherical sym- 



metry, we approximate the host profile by averaging its 
density in spherical shells around the satellite location. 
See Fig. [5| Thus, for the average host profile we have 
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where = + — 2drcos9, d is the center-to-center 
distance between the halos, and 9 is the azimuthal angle. 
Averaging the satellite density psat around its own center 
leaves its own NFW profile unchanged. We insert the 
total density 



P = Psat + (Ph) 



(21) 



in Eq. (14) with the boundary conditions (15), yielding 



the radial symmetron profile. The Newtonian force is 



given by integrating the profile Eq. ( 20 ) in spherical shells 



to find the enclosed mass (see Appendix \M. We note 



that this averaged density profile will slightly overesti- 
mate (underestimate) the force deviation F^/F^ on the 
side of the satellite nearer to (farther from) the host. 

Figure |6] shows the density profile of a cluster-size (2 x 
W^Mq) host modeled by Eq. ([20|, an NFW sateUite, 
and their sum, for two host-satellite separation d — 1 and 
2 .Rhost- In the inner parts of the satellite its own NFW 
profile dominates the total density, while the host profile 
is a slowly increasing function of r. At r ~ i?sat the 
averaged host density becomes the dominant component. 
After the cusp, which occurs at the host center we see the 
host density rapidly transitions to the NFW p cx 
power law. The pictured density profile is qualitatively 
similar for all relevant satellite and host masses. 



B. Results 



The cluster's screening effect on the satellite is evident 
in Fig. T\ At low radii the force deviation F(f,/F^ for 
the totardensity profile has a similar slope to that of the 
isolated satellite profile, but with an amplitude decreased 
below the observationally interesting level of ~ 1%. In 
each panel there is a sharp transition of the solid line at 
the distance d (i.e., location of the cusp in Fig. [6]) sepa- 
rating the host and satellite. Here the force ratio for the 
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FIG. 7: Force deviation F^/Fisi for the total density profile of Eq. (21 1 (solid lines). This is the deviation that would be 
experienced by a test mass in the neighborhood of a satellite halo which is itself blanket screened by a nearby host (A/host = 
2 X IQ^'^Mq). The deviation is shown for two values of the host-satellite separation d: thick solid lines for d = 1-Rhost and thin 
solid lines for d = 27?host • Also shown are / Fn for isolated halos with mass equal to the satellite (dotted) and host (thin and 
thick dot-dashed). The environmental screening from the host brings the modified forces of the satellite below 10%, therefore 
nearly unobservable in each case. Figure |8] shows larger deviations for lower host masses. 



total density profile is closely approximated by that of 
the host alone. At very large radii, all three curves fall 
off rapidly, as the symmetron force decays exponentially 
due to its finite compton wavelength. The environmen- 
tal screening from the host brings the force deviation of 
the satellite below the level of observations for each case. 
Note that for this and the following two-halo calculations, 
we set (f)o = M^/Mpi. 

In order to see the partial screening of a satellite, it 
is necessary to consider smaller hosts. Fig. [8] shows 
the force deviation F^/F^ for the same range of satel- 
lite masses and host satellite separations as in Fig. [g but 
with host masses of 10^^ and lO^^M©. For Msat < Mhost, 
including the host at 2 i?host decreases the fifth force by 
50-150% relative to gravity, and halving the separation to 
1 i?host brings a further decrease of 5-10%. For example, 
an isolated lO^^M© has sa 150%Fn at its virial radius, 
while if a lO^^M© halo is 2 or 1 virial radii away, this de- 
viation is cut down to 60%Fn or 45%^n, respectively. 
See Appendix |B] for calculations of the modified forces at 
the satellite virial radius as a continuous function of d. 

A few comments on the validity of the approximate 
host profile of Eq. (20) and the results of Figs. [Tllslare 



in order here. Since the NFW profile of the host does not 
change much across the diameter of the smaller satellite, 
we expect the force deviation predictions to be accurate 
in the regime r < i?sat- Here our spherically symmet- 
ric approximation should slightly overestimate the devi- 
ation on the side of the satellite nearer the host, while 
underestimating the deviation on the far side. For r ^ d 
the offset between host and satellite is negligible, so that 
Eq. (20) approaches the NFW density profile Eq. (11). 



In the intermediate regime r ~ d we note that the 
screening will vary widely between the near and far sides 
of the satellite. On the side nearer the host we should 
observe thorough blanket screening, while on the far side 
screening from the host may well be negligible, depend- 
ing on the separation d. However, given that dynamical 
tracers of the satellite mass (such as stars, HI gas, and 
satellite galaxies) are generally confined to r < i?sat: our 
model describes well the observationally relevant range. 



V. DISCUSSION 

We have presented predictions for the modified forces 
of symmetron gravity around dark matter halos, model- 
ing the halo density with the NFW profile. For a large 
range of halo masses 10""^° — 10^^ Mq we find order unity 
deviations from GR at distances 1 — 7 i?3oo from the 
halo center, while the dwarf-size halos of 10^° — IO^^Mq 
exhibit large deviations throughout (Fig. [2]). 

We have also shown an average of this deviation given 
by the quantity ^i^^.^. suggested by P5; our theoretical cal- 
culations of ^^.^ (Fig. 3| are simply related to observables 
by Eqs. ^ and 0f We find > 20% differences over 
10^^ Mp), indicating that observa- 



Thus, our approximation captures well the small and 
large r behavior of the exact host-satellite system. 



the mass range 10 — iwiviq, 
tions of dynamical and lensing masses of galaxies are a 
promising way to constrain parameters of the theory. 

Furthermore, ^^j^ is necessary to test the mass function 
of the theory. Reference |19| obtains the mass function 
from symmetron N-body simulations by counting halos 
as a function of their true or lensing mass. However, 
observations more commonly yield the dynamical mass of 
virialized structures; a means of converting between the 
two is therefore essential to constrain the theory using 
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FIG. 8: Force deviation F^/Ft^ for the total density profile of Eq. (21 1 (solid lines), using smaller host masses than in Fig. 
This is the deviation that would be experienced by a test mass in the neighborhood of a satellite halo which is itself (partially 
blanket screened by a nearby host. The deviation is shown for two values of the host-satellite separation d: thick solid lines 
for d = 1-Rhost and thin solid lines for d = 27?host. For comparison, the dotted line plots the deviation for an isolated satellite: 
the difference between dotted and solid lines indicate the additional screening due to the nearby host. In each case, the host 
significantly decreases the force deviations of the satellite, but to a level that is still potentially observable. 



the mass function. 

We have gathered predictions for symmetron, 
chameleon, and Vainshtein screening in the neighbor- 
hood of NFW haloes in order to find ways to distinguish 
between classes of MG models in realistic astrophysical 
situations. We find significant differences among the 
three screening mechanisms, including 

• In contrast to the chameleon, the lack of a dis- 
tinct thin-shell radius for the symmetron results 
in a nonzero deviation at small radii even for par- 
tially screened halos. Since visible tracers of galax- 
ies are often well within the virial radius, the tests 
of [211 [221 125] are more easily applied to test sym- 
metron screening. 

• At a given radius, the chameleon deviations change 
rapidly with mass and therefore reach their maxi- 
mum of 4/3 quickly. The symmetron exhibits more 
gradual changes in the deviation, while DGP mod- 
els have no dependence on halo mass. 



We note that there is an approximation involved in us- 
ing the NFW profile: although this profile is a good fit 
to stacked halo densities in GR simulations, it is possi- 
ble that the profile itself will change with the modified 
forces. However, the symmetron N-body simulations of 
[l9] are able to resolve halos of mass > 5 x lO^^M^g; for 
these larger halos we have compared our results for the 
symmetron field and find them qualitatively consistent 
with the N-body simulations. (Note that [H] defines a 
parameter zssb related to ji and Ms (see [19] for details) 
and focus on 2;ssb = 2.0. We have checked our results 
with theirs using this same value for zssb, but our stan- 
dard choice of parameters corresponds to zssb ~ 0.5, 
for which the theory predicts relatively lesser deviations 
from GR. Thus, consistency with the simulations for the 
larger value of zssb is actually a more stringent test than 
consistency for our choice of parameters.) 

This work has not considered two other screening 
mechanisms that have recently been discovered: the en- 
vironmentally dependent dilaton [SHI HZ] and k-mouflage 
[25] . However, k-mouflage screening has been shown to 
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be similar to Vainshtein (i.e., independent of halo mass), 
and due to the similar dependence of the derivative of the 
coupling {dA/d(j) oc (j)) in symmetron and dilaton mod- 
els, the dilaton screening may exhibit similarities to the 
symmetron. 

We have also shown results for screening of satellite ha- 
los by larger neighboring hosts. We have approximated 
the host density by averaging its NFW profile in spherical 
shells about the satellite: the resulting profile is shown 
in Eq. ( 20 1 and Fig. |6] We found that the environmental 
screening effect from a 2 x 10^"* cluster is sufficient to 
reduce deviations from GR well below the level of 1% for 
smaller halos located within twice the host virial radius 
(Fig. [t]). However, Fig. |8] shows that for smaller hosts 
of mass 10^^ — the force deviations around the 

satellite are decreased (relative to the isolated case) but 
still potentially observable. One caveat to this method 
should be mentioned: in reality, adding a second halo 
breaks the spherical symmetry of our system. Thus we 
expect that our results for two-body systems slightly un- 
derestimate (overestimate) the screening on the side of 
the halo nearer to (farther from) the second halo. 

We comment briefly on the possibility of analytical 
solutions for two-body systems in symmetron modified 
gravity. There has been some success in finding ap- 
proximate analytical solutions of two-body systems in 
chameleon theories. In the thin shell effect is used 
to find solutions for uniform density spheres in a variety 
of configurations. Furthermore, 30 expands these solu- 
tions to include an NFW host halo and test body (point 



mass) satellites of Mgat < \Q^^Mq. However, these re- 
sults do not straightforwardly translate to symmetron 
gravity: the thin-shell effect of the chameleon allows a 
clean division into 3 regions in the neighborhood of an 
object with a thin shell (inside the shell, the thin shell 
itself, and outside the shell); in contrast, the symmetron 
transitions more smoothly within screened objects (see 
Fig. |4]), making such a partitioning much more difficult. 
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Appendix A: Newtonian force of host 

Here we show the result of integrating the density pro- 
file of Eq. ( 20 1 to find the Newtonian force of the host 
component. We have 
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Appendix B: Variations of host-satellite separation 



Although using the averaged host profile of Eq. (20) 



preserves spherical symmetry, we must still solve the non- 
linear Eq. ( 14 1 for each choice of host-satellite separation 
d. To see the effect of varying this parameter continu- 
ously, we make a different approximation. If the sym- 
metron profile of the host at the location of the satel- 
lite has value 0host(c?), and this value varies little across 



the diameter of the satellite, we can instead solve the 
approximately equivalent system of an isolated satellite 
with asymptotic field value 4>{r — c») = (t>host{d). 

Fig. [9] plots the satellite force deviation at li?sat as a 
function of d. This is done using two different approxi- 
mations: the method using the asymptotic value of the 
field appropriate for a small object at distance d from 
the host and the total density profile method of Eq. ( 20 1 . 



Plotting continuous curves for the latter is impractical: 
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FIG. 9: (left panel): Satellite force deviation F^/F]ss at 1 i?aat as a function of center-to-center distance d from a nearby host halo. 
Continuous lines display the result using the asymptotic value of the field appropriate for a relatively small object at distance d 
from the host (see text for details). Satellite masses Msat ~ 10^^, 10^^, W^^Mq correspond to solid, dashed, and dotted curves, 
respectively. Also shown for the same range of masses are the results using the total density profile of Eq. ( |20[ ) applied at a 
discrete set of points. Horizontal lines show the force deviation of isolated satellites. The host mass is 2 x 10 A/q. For the 
10^^ and 10^'^Mq satellites, the screening from the host has a significant (~ 10%) effect even at separations of ~ 7 — 8 -Rhost. 
(right panel): Same as left, but for a host mass of IO^^Mq. 



the shooting method must be employed to solve Eq. ( 14 ) 



for each distance d. Fortunately, a discrete set of points 
shows the trend sufficiently well. 

We see that the computationally much simpler asymp- 
totic value method approximates well the more sophisti- 
cated average density calculation as long as the host and 



satellite masses are within 2 orders of magnitude. Here 
the simpler method overestimates the fifth force by only 
5 — 30% Fn over the range of separations considered. The 
differences exceed this level for the case of a IO^^Mq halo 
in the neighborhood of a 2 x 10^'* host, but only at 
separations of 5-6 i?host where it reaches « 50%. 
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